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Abstract 

The MSW (Mikheyev-Smirnov-Wolfenstein) effect is the adiabatic or partiahy adi- 
abatic neutrino flavor conversion in medium with varying density. The main notions 
related to the effect, its dynamics and physical picture are reviewed. The large mix- 
ing MSW effect is realized inside the Sun providing the solution of the solar neutrino 
problem. The small mixing MSW effect driven by the 1-3 mixing can be realized for 
the supernova (SN) neutrinos. Inside the collapsing stars new elements of the MSW 
dynamics may show up: the non-oscillatory transition, non-adiabatic conversion, time 
dependent adiabaticity violation induced by shock waves. Effects of the resonance 
enhancement and the parametric enhancement of oscillations can be realized for the 
atmospheric and accelerator neutrinos in the Earth. Precise results for neutrino oscil- 
lations in the low density medium with arbitrary density profile are presented and the 
attenuation effect is described. The area of applications is the solar and SN neutrinos 
inside the Earth, and the results are crucial for the neutrino oscillation tomography. 



*Talk given at the Nobel Symposium 129, "Neutrino Physics", Haga Slott, August 19 - 24, 2004. 



1 Introduction 



Variety of matter effects [1] on propagation of mixed neutrinos [2, 3] depends on (i) channel 
of mixing - the type of neutrino states involved: active, sterile, mass eigenstates, etc.; (ii) 
density profile: constant, monotonous, periodic, fiuctuating; (iii) properties of medium: 
polarization, motion, chemical composition (thermal bath and neutrino gases are special 
cases); (iv) presence of new non-standard neutrino interactions. 

Dynamics of the matter effects can be classified by degrees of freedom involved. In the 
2z/ case, an arbitrary neutrino state can be expressed in terms of the eigenstates of the 
instantaneous Hamiltonian, and 1/2^1, as 

u{x) = cos 9al^lm + Siu Oaiy2me'^"' , (l) 

where 

• &a — dai^) determines the admixtures of eigenstates in ^{x); 

• ^rn{x) is the phase difference between the two eigenstates (phase of oscillations): 

<^m{x)= r AHdt', (2) 

Jxo 

here AH = — Him is the difference of eigenvalues. Eq. (2) gives the adiabatic 
phase, an additional contributions can be related to non-adiabaticity and topology. 

• The mixing angle in matter, defined as {ve\^im) = cos^^, {i'e\^2m) = sin^^, etc., 
determines the flavor contents (or flavors) of the eigenstates. The angle is the function 
of matter density ne{x): 9m — 9m{ne{x)). 

Besides these, the effects of absorption and loss of the coherence can change the nor- 
malization of the eigenstates or suppress interference. 

Different processes are associated with these degrees of freedom. In particular, in pure 
form the oscillations [2, 3, 4, 1, 5] are the effect of the monotonous phase difference increase 
$m, which occurs in the uniform medium. In contrast, the MSW effect [1, 6] is associated 
to the change of flavors of the eigenstates (change of 6m{x)) in the nonuniform medium. In 
general, an interplay of different effects occurs which is induced by simultaneous operation 
of several degrees of freedom. 

In this paper I will discuss only few effects selected on the following ground: usual 
massive neutrinos with the flavor mixing and standard weak interactions; the effects in the 
Sun, collapsing stars and the Earth; the effects which are reahzed in Nature or have a good 
chance to be realized and established in our future studies. 
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2 The MSW effect. 



The MSW effect [1, 6] is the adiabatic or partially adiabatic neutrino flavor conversion in 
medium with varying density. The flavor of neutrino state follows the density change. Here 
I review the main notions related to the effect and describe its dynamics. 



2.1 Main notions 



1. Refraction. At low energies the elastic forward scattering only is relevant in most of 
applications [1, 7]. It can be described by the potential, V^, or equivalently, the refraction 
index: n^e/ — 1 = V /p. 

The difference of potentials (for a single neutrino) in different spatial points V = V{x) 
leads to the "neutrino optics" phenomena such as reffection, complete internal reflection, 
banding of neutrino trajectories, focusing by the stars, etc.. The values of refraction index 
are very close to unit, e.g..,'. n^ef — 

1 = 10-2° inside the Earth, 10"^^ inside the Sun, lO^^ 
in the neutron stars. So, the "neutrino lenses" should be of the astrophysical size. 

The difference of potentials for different neutrinos, z/g and u^. V = Ve — Va influences 
evolution of mixed neutrinos (even for constant potentials). In usual matter the difference 
is due to the charged current scattering of z/g on electrons {u^e — >• u^^) [1]: 

V = Ve-Va^ V2GFn, , (3) 

where G^? is the Fermi coupling constant and is the number density of electrons. It 
leads to appearance of additional phase in the neutrino system: A(j)matter = (K — Va)t- The 
distance over which this "matter" phase equals 27r determines the refraction length: 

2. Eigenstates and mixing in matter. In the presence of matter the Hamiltonian of system 
changes: Hq ^ H = Hq + V, where Hq is the Hamiltonian in vacuum. Correspondingly, 
the eigenstates and the eigenvalues of H change: z/i,z/2 '^im, ^2m, m\/2E,m\/2E — > 
Him-i H2m- Here ui and U2 are the mass eigenstates with masses m\ and m2. 

The mixing in matter is defined with respect to the eigenstates and V2m- Similarly 
to the vacuum case, the mixing angle in matter, 9m-, determines relations between the 
eigenstates in matter and the flavor states: = cos 9mi^im + sm9mt^2m, etc.. In matter, 
both the eigenstates and eigenvalues, and consequently, the mixing angle depend on matter 
density and neutrino energy. 

3. Resonance. Dependence of the effective mixing parameter in matter, R = sin^ '^9^, 
on ratio of the vacuum oscillation length, Zj, = AirE/Am'^, and the refraction length: x = 
lu/lo = 2EV/ Arn^ oc Ene (fig. 1) has a resonance character [6]. At 

lu — h cos 29 (resonance condition) (5) 
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Figure 1: The dependence of the effective mixing parameter R — sin ori x — Iu/^q Eti^ for 
two different values of the vacuum mixing: sin^ 26 = 0.825 (black, red) which is the LMA mixing 
and tan^ 6 = 0.08 ( grey, green). The semi-plane x < corresponds to the antineutrino channel. 

the mixing becomes maximal: sin^ 26m = 1- For small 6 the condition (5) reads: 

Vacuum oscillation length ^ Refraction length. (6) 

That is, the eigenfrequency which characterizes a system of mixed neutrinos, 1/lu, coincides 
with the eigenfrequency of medium, I/Iq. For large vacuum mixing (LMA has cos 26' = 0.4) 
there is a significant deviation from the equality (6). Large mixing corresponds to the case 
of strongly coupled system for which the shift of frequencies occurs. 

The resonance condition (5) determines the resonance density: 

= no cos 2^, no = iV 

The width of resonance on the half of height (in the density scale) is given by 

2 Anf = 2nf tan 26 = 2no sin 26. (8) 

When the vacuum mixing approaches maximal value, 6 = it/A, the resonance shifts to zero 
density, — >■ 0, and the width of resonance Anf increases converging to ng. 

In medium with varying density, the resonance layer is determined by the interval in 
which the density changes from n^ — Anf to + An^. 

4- Adiabaticity. Since in non-uniform medium the density changes on the way of neutrinos, 
ne = ne{x), the Hamiltonian of system depends on time (distance): H = H{t). Therefore, 
(i) the mixing angle changes in course of propagation: 6m = 6m{ne{x)); (ii) the eigenstates 
of instantaneous Hamiltonian, uim and z/2m, are no more the "eigenstates" of propagation, 
and the transitions uim ^ i'2m occur. 
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If the density changes slowly, the system (mixed neutrinos) has time to adjust the change 
leading to the adiabatic evolution [1, 6, 8, 9]. The adiabaticity condition is [9] 
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H2m — H 



< 1. 



(9) 



Im 



As follows from the evolution equation for the neutrino eigenstates [6, 9], \6„i\ determines 
the energy of transition ^ v^mi and \H2m — -f^iml gives the energy gap between the 
levels. So, the condition (9) means that the transitions uirn ^ ^2m can be neglected and 
the eigenstates propagate independently {6a — const in (1)). 

If 9 is small, the adiabaticity is critical in the resonance. It takes the form [6] 



Arfl > Ir, 



(10) 



where Ir — lyj sin 2^ is the oscillation length in resonance, and IStr — nftaiXi29/{dne/dr)R 
is the spatial width of resonance layer. The adiabaticity condition has been considered out- 
side the resonance and in the non-resonance channel in [10]. In the case of large vacuum 
mixing the point of maximal adiabaticity violation, n"*^ [11], is shifted to densities larger 
than the resonance one: n'^ no > n^. 



2.2 Dynamics of the MSW effect. 

1. Dynamical features of the effect in the adiabatic case can be summarized as follows. 

• The flavors of eigenstates change according to density change; the flavors are deter- 
mined by 9m{x) = 9m{ne{x)). 

• The admixtures of the eigenstates in a propagating neutrino state do not change due 
to the adiabaticity; there is no Vim ^ t'2m transitions. The admixtures are given by 
the mixing in production point, 9^. 

• The phase difference ^rn{x) (2) increases leading to oscillations. 



Two degrees of freedom are operative: the phase <I>m and the flavor, 9„i. The MSW effect 
is driven by the change of flavors of the neutrino eigenstates in matter with varying density. 
The change of phase produces the oscillation effect on top of the adiabatic conversion. 

2. Spatial picture of the MSW effect. In fig. 2 shown are dependences of the average 
probability, P, and depth of oscillations given by P"*"^ — p"^*"^ on 
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Figure 2: Dependences of the average probability (dashed hne) and the depth of oscillations 
^pmax ^ pmm _ gQjj(j lines) on y for yo = ~5. The resonance layer corresponds to y = and the 
resonance layer is given by the interval y = —1-^1. For tan^ 6 = 0.4 (large mixing MSW solution) 
the evolution stops at yj = 0.47. 

the distance (in the density scale) from the resonance density in the units of the width of 
resonance layer [6]. In terms of n the conversion pattern depends on the initial value uq 
only which reflects universality of the adiabatic evolution. The probability is the oscillatory 
function which is inscribed into the band shown by the solid lines. There is no explicit 
dependence on the vacuum angle 9. With decrease of yo, the oscillation band becomes 
narrower approaching the line of non- oscillatory conversion. For zero final density: yj = 
1/ tan 2^. The smaller the mixing (and therefore, the larger final ?//) the stronger transition. 

3. Adiabaticity violation. If density changes rapidly, so that the condition (9) is not 
satisfied, the transitions ^ i^2m become efficient [6, 12]. Therefore admixtures of the 
eigenstates in a given propagating state change: 9a = 9a{x). Now all three degrees of 
freedom - phase, flavor, and admixture - become operative. Typically, adiabaticity breaking 
leads to weakening of the flavor transition and enhancement of oscillations. 

4. Graphic representation [13] is based on the analogy of the neutrino evolution with the 
behavior of spin of the electron in the magnetic fleld. The neutrino evolution equation can 
be written as 

f.(Bx.). (12, 
where the neutrino vector of length 1/2 (equivalent of spin) is 

u={R,I,P- 1/2) = (Rei/lz/^, Imulu^, v\v, - I/2) , (13) 

'^a = i^a{x), (a = e) are the neutrino wave functions, and P is the survival probability ^. 
^The elements of this vector are nothing but components of the density matrix. 
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Figure 3: Graphic representation of the MSW effect. The direction of the cone axis is determined 
by 26„i, the cone angle is given by 6a and the position of the neutrino vector on the surface of 
the cone is fixed by the phase In the adiabatic case (left panel) the direction of the axis 

flips according to ^rn^^e change, and the cone angle is unchanged. In the non-adiabatic case (right 
panel) also the cone angle changes. 

The vector of "magnetic field" equals 



where Im — 27r/Aif is the oscillation length in matter. The vector v moves on the surface 

— * 

of the cone with axis B according to increase of the oscillation phase, 



General conditions for the MSW conversion are: (i) slow enough density change; (ii) crossing 
the resonance layer; (iii) large enough matter width (minimal width condition) [14] . These 
conditions are satisfied for the solar neutrinos inside the Sun, for supernova neutrinos inside 
collapsing stars and can be satisfied for neutrinos in the Early Universe. 

3.1 Solar Neutrinos. Large Angle MSW solution 

The large mixing MSW conversion provides the solution of the solar neutrino problem [15] . 
The best fit values of the oscillation parameters from combined analysis of the solar and 
KamLAND data (in assumption of the CPT invariance) are [16]: 




2-K 



(14) 



3 Realizations of the MSW effect 



= 7.9 . 10-5 eV^ tan^ 9 = 0.40. 



(15) 
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Analysis of the solar neutrino data alone leads to smaller mass split, Am^ = 6.3 • 10 ^eV^ 
which agrees with (15) within la. 

1. Physical picture. According to LMA, inside the Sun the initially produced electron 
neutrinos undergo the highly adiabatic conversion: z/g — > cos^^^^z^i + s\n6^^V2^ where 9^ is 
the mixing angle in the production point. On the way from the central parts of the Sun 
the coherence of neutrino state is lost after several hundreds oscillation lengths [17], and 
incoherent fluxes of the mass states i/i and V2 arrive at the surface of the Earth. In the 
matter of the Earth vi and z/2 oscillate partially regenerating the i/e-flux. The averaged 
survival probability can be written as 

P = sin^ e + cos^ cos 29^2 - cos ^O'^^Ueg, (16) 

where the first term corresponds to the non-oscillatory transition (dominates at the high 
energies), the second term is the contribution from the averaged oscillations which increases 
with decrease of energy, and the third term is the regeneration effect /^gg. At low energies 
P reduces to the vacuum oscillation probability with very small matter corrections. 

2. Status of LMA. The solution provides a very good global fit of the solar neutrino data. 
There is no statistically significant deviation from description given by the standard solar 
model (SSM) [18] and the LMA solution. 

The key observation which testifies for the MSW (matter) effect in the Sun is stronger 
than 1/2 suppression of the signals at SK and SNO. The fg-survival probability extracted 
from the CC/NC ratio at SNO is 

(P,e) = 0.3lt°:^^ (3a), (17) 

that is, (Pee) < 0.50 at 5(7, whereas the vacuum 2z/ oscillations can produce (Pee) > 0.5. 

Observations of two other signatures of the solution (i) the upturn of the spectrum at 
low energies {E < 7 — 8 MeV), (ii) the day-night asymmetry of signals with larger fiux 
during the night are the main objectives of the forthcoming studies. 

No viable alternative to the LMA solution exists and possible effects beyond LMA are 
substantially restricted already now. 

There is a very good agreement of the results from solar neutrinos and KamLAND 
which imphes the CPT conservation. Furthermore, it shows correctness of theory of both 
the vacuum oscillations and conversion in matter. 

3. Testing the MSW effect in the Sun. Important way to test the effect is to introduce 
the free parameter, aMSW, in the the matter potential 

V ^ auswV, (18) 

and to determine gmsw from the data [19]. The global analysis of the solar and reac- 
tor (KamLAND -|- CHOOZ) results with Aml^, sin^^i2 and ausw being unconstrained 
gives [19] ausw — 0.8 — 3.0 (95% C.L.) with the best fit value ausw — 1-6; zero value of 
dMSw is excluded at 6a. 
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4- Precision measurements in solar neutrinos. Identification of the LMA solution opens 
new possibilities in [20] (i) precise description of the LMA conversion both inside the Sun 
and in the Earth taking into account various corrections; (ii) estimation of accuracy of 
approximation made; (iii) obtaining the accurate analytic expressions for probabilities and 
observables as functions of the oscillation parameters. There are three small quantities 
which allow for a very precise expansions. 

(1) . Smallness of the adiabaticity parameter 

where h is the height of solar density profile, allows to use the adiabatic perturbation theory. 
The non-adiabatic corrections to the averaged survival probability are of the order 7^ [20]. 

(2) . Smallness of the ratio Arprod/h, where Ar^rod is the size of the neutrino produc- 
tion region, allows one to make the averaging of the survival probability over the neutrino 
production region in the analytic form [20]. 

(3) . Smallness of the parameter 

2EVf 

where Ve is the matter potential in the Earth, allows one to develop a very precise pertur- 
bation theory for the neutrino oscillations inside the Earth [20]. 



3.2 MSW effect and Supernova neutrinos 

In supernovae one expects new elements of the MSW dynamics. The SN neutrinos probe 
whole level crossing scheme, and the effects of both resonances (due to Am^2 ^"^13) 
should show up. Various effects associated to the 1-3 mixing can be realized, depending 
on value of ^^13 (fig. 4). As follows from fig. 4, the SN neutrinos arc sensitive to sin^ ^^13 as 
small as 10~^. Studies of the SN neutrinos will also give information on the type of mass 
hierarchy [21, 22, 23, 24, 25]. 

1. The small mixing MSW conversion. This can be realized due to the 1-3 mixing and 
the "atmospheric" mass split Amfg. 

2. The non-oscillatory adiabatic conversion [6] is expected for sin^ 9is > 10~^. The density 
in the production point is extremely large: n(0) ^ n^, and therefore the mixing in the 
initial state is strongly suppressed. So, the neutrino state coincides with the cigcnstatc in 
matter: z/g = i'2m- Due to adiabaticity (no z/j^ ^ i^jm transition) the neutrino state coincides 
with this eigenstate during whole evolution: i/(t) = i/2m- The interference is negligible since 
simply there is no second component to interfere with, and consequently, oscillations are 
absent. The flavor of the neutrino state changes as the flavor of the eigenstate and the 
latter follows the density change. 
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Figure 4: Scales of the 1-3 mixing probed with supernova neutrinos. Indicated are the regions of 
various effects in neutrino propagation. 

3. Adiabaticity violation occurs if the 1-3 mixing is small sin^ Oi^ < 10^^. 

The shock wave can reach the region of the neutrino conversion, p ~ 10^ g/cc, after 
= (3 — 5) s from the bounce (beginning of the u— burst) [26]. Changing suddenly the 
density profile and therefore breaking the adiabaticity, the shock wave front influences the 
conversion in the resonance characterized by Amfg and sin^^ia, if sin^^is > 10"'^. 

The following shock wave effects can be seen in neutrinos (antineutrinos) for normal 
(inverted) hierarchy: (1) change of the total number of events in time [26]; (2) wave of 
softening of the spectrum which propagates in the energy scale from low energies to high 
energies [27]; (3) delayed Earth matter effect in the "wrong" channel {e.g., in neutrino 
channel for normal mass hierarchy) [21]. Modification of the density profile by the shock 
wave leads to appearance of additional resonances below the front [28]. Reverse shock 
produces a "double dip" time feature in the average neutrino energy [29, 25]. Monitoring 
the shock wave with neutrinos can shed some light on the mechanism of explosion. 

4. Neutrinos from SN1987A. After confirmation of the LMA MSW solution we can defi- 
nitely say that some effect of flavor conversion has already been observed in 1987. 

In the case of normal mass hierarchy the adiabatic z/g — > i^i and Ufj,^r —>■ ^2 transitions 
occurred inside the star, and then ui and z/2 oscillated inside the Earth [10, 21]. In terms 
of the original fluxes of the electron, and muon antineutrinos, F^{ue) and F°(z/^), the Ue— 
flux at the detector can be written as 

F(z7e) = F0(z7e)+pAF°, (21) 

where AF^ = F{i'^) — F{ue), and p is the permutation factor which can be calculated 
precisely: p = sin^ — freg- Due to difference in the distances traveled by neutrinos to 
Kamiokande, 1MB and Baksan detectors inside the Earth, the Ue regeneration factors freg 
and therefore p differ for these detectors. This can partially explain the difference of the 
Kamiokande and 1MB energy spectra of events [30]. 

One must take into account the conversion effects in analysis of SN1987A [30] as well as 
future supernova neutrino data. The conversion can lead to increase of the average energy 
of the observed events by (30 - 40)%. Inversely, not taking into account the conversion effect 
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produces errors in determination of the average energy of the original Ve spectrum up to 40 
- 50 % in Kamiokande, and factor of 2 in 1MB. 

For the inverted mass hierarchy and sin^ > 10^^ one would get nearly complete 
permutation, p 1, and therefore a harder z/g spectrum, as well as an absence of the Earth 
matter effect. This is disfavored by the data [10, 22], see however [23]. 

4 Matter effects in Neutrino Oscillations 

Pure oscillation effect can be realized in the uniform medium. Mixing is constant, 6m{E, n) = 
const., and therefore 

• the flavors of eigenstates do not change; 

• the admixtures of eigenstates do not change; there is no i^i^ transitions: i^i^ 
and are the eigenstates of propagation; 

• monotonous increase of - the phase difference between the eigenstates occurs. 

Only one degree of freedom operates - the phase $,„ and all others are frozen. This is 
similar to what happens in vacuum. The depht and length of oscillations are determined 
by the mixing and energy splitting in matter: sin^ 2^j„, 1^ — 2t:/ {H2m — Hi^). 

The oscillations are reahzed in the Earth which can be considered as the multi-layer 
medium with nearly constant density in each layer. Variety of possibilities exists depending 
on the neutrino trajectory (zenith angle), neutrino energy and channel of oscillations. 

4.1 Resonance enhancement of oscillations 

For a given (constant) density, the length and depth of oscillations depend on the neutrino 
energy. This leads to a characteristic distortion of the energy spectrum, F{E)/ Fq{E), where 
Fo(E) and F{E) are, e.g. the spectra of {e.g. v^) neutrinos in the source and detector 
correspondingly (flg. 5). The ratio F{E)/Fq{E) given by the survival probability has an 
oscillatory dependence on the energy. At the resonance energy, Er, the oscillations proceed 
with maximal depth; they are enhanced in the resonance range [6]: E = Er ± AEr, where 
AEr = tan2eER = sin 2^^^ and = AmV2V2G'i.ne. 

The effect is realized for the high energy neutrinos {E ~ GeV) in the mantle of the 
Earth: constant density is a good flrst approximation. The effect is expected to be seen 
in the accelerator long baseline experiments [31], and future high statistics atmospheric 
neutrino studies [32, 33]. 
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Figure 5: Resonance enhancement of oscillations in matter with constant density. Shown is the 
dependence of ratio of the final and original fluxes, (of e.g. Ve) F/Fo, on energy (x = lu/lo oc E) for 
a thin layer, L = Iq/tt, and sin^ 26 = 0.824 (left panel), and thick layer, L = lOlo/ir, and tan'^ 6 = 
0.08 (right panel). The oscillatory curves are inscribed in to the resonance curve (1 — sin'^29„i). 



4.2 Parametric enhancement of oscillations 

This enhancement is related to certain condition for the phase of oscillations [34, 35]. It 
provides another way of getting strong transition: no matter enhancement or resonance 
conversion are needed. No large or maximal mixings in vacuum or matter are required. 

The simplest case which can be realized in Nature is neutrinos in the castle wall profile. 
The latter consists of the alternate layers with two different densities [35, 36, 37]. Let $i, 
$2 and 6"^, 6^ be the phases and mixing angles in the layers 1 and 2. Then under condition: 

SiC2Cos2^™ + S2CiCos2^^ = 0, (22) 

where Si = sin$j/2, ci = cos$j/2, {i = 1,2), which is called the parametric resonance 
condition [38], the flavor transition can be complete. Simple realization of (22) is 

$1 = $2 = TT (23) 

which leads to Ci = C2 = 0. Eq. (22) can be satisfied for neutrinos (the z/^ — Uf. channel 
with the Am^3 and 1-3 mixing) with few GeV energies which cross the core of the Earth 
[32, 37]. These neutrinos propagate in three layers of matter: mantle-core-mantle. In the 
approximation of constant densities inside the layers, the profile can be considered as a part 
of the castle wall profile [36]. 

For small mixings in both layers, the maximal enhancement occurs when the condition 
(23) is satisfied (fig. 6, left panel). Apparently for few more layers that would lead to maximal 
transition. On the other hand, even three layers are enough to get nearly maximal transition 
provided that the mixings in matter are not small (fig. 6, right panel). 

For Am^3 = 2 ■ 10"'^ eV^, the strongest parametric effect (fig. 6 right panel) can be 
observed in a sample of the atmospheric neutrinos with E > 2.5 MeV (2 times larger than 
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Figure 6: Graphic representation of the parametric enhancement of neutrino oscihations. The 
notations are explained in fig. 3 and in Eqs. (13, 14). The left panel: the neutrino energy above 
the resonance energy in the mantle (matter suppressed mixing). The right panel: the neutrino 
energy is between the core and mantle resonance energies. The mixing in the core is large; in the 
mantle $1 > tt/A and in the core > 27r. 

the energies of multi-GeV sample [32]). Manifestation is the excess of the e-like events for 
the core crossing trajectories. It can be seen as an enhanced up-down asymmetry of the 
e-like events [32, 33]. 

4.3 Neutrino oscillations in the low density medium 

The condition of low density. 



means that the potential energy is much smaller than the kinetic energy. In this case one 
can use small parameter e(x) (20) to develop the perturbation theory [39]. 

For the LMA oscillation parameters and the solar and supernova neutrinos: e(x) = 
(1 — 3) ■ 10~^. The relevant channel of oscillations is the mass-to-flavor transition, V2 ^ei 
since both the solar and SN neutrinos arrive at the Earth as the incoherent fluxes of mass 
eigenstates. The probability of V2 — > v^. can be written as = sin^ 6' + /reg, and oscillations 
appear in the first order in e. Using the e— perturbation theory the following expression for 
the regeneration factor f^eg has been obtained [39, 40] 



Here Xq and xj are the initial and final points of propagation correspondingly, $m(a; ^ ^f) 
is the adiabatic phase (2) acquired between a given point of trajectory, x, and final point, 
x/. The latter feature has important consequence leading to the attenuation effect. 




(24) 




(25) 
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In the second order in V, and therefore e, the regeneration factor becomes [41] 



freg = Slli^ 29 



sin $r 



[Xr 



]I + cos29 f 



(26) 



where for the symmetric profile 

I = 



dxV{x) cos 



[Xr 



x). 



(27) 



In / the integration proceeds from the central point of trajectory, Xc, to the final point, Xf. 
The phase is calculated from Xc to a given point x. Essentially, the integral / plays the role 
of expansion parameter and it can be estimated as / < 2EVmax/^m'^ = ^max- 

The perturbation theory [41] can be improved if the expansion is performed with respect 
to certain potential Vq rather than zero. The effective expansion parameter becomes smaller. 



4.4 Attenuation effect 

Eq. (25) allows one to estimate sensitivity of the oscillation effects to structures of the 
density profile [39]. Consider some structure in the point x of the trajectory at the distance 
d = Xf — X from the detector. According to (25) for the mass-to-fiavor transition the 
potential V{x) is integrated with sin $„((/). The larger d, (and therefore, the larger $m((i)), 
the stronger averaging effect when some integration over the energy is performed. So, weaker 
sensitivity of the oscillation effects to remote structures of the profile should show up. The 
integration over energy with the energy resolution function R{E, E') , 

freg{E) = j dE' R{E , E') Ug{E') , (28) 

can be expressed as 

- 1 f'^f 

freg^ 7i^^^'^29 dxV{x)F{xf-x)sm^rn{x^Xf), (29) 

Jxn 



where F{xf — x) is called the attenuation factor [39]. For the box-like function R{E,E') 
with width (energy resolution) AE we obtain 

^/ ,x 1 . TidAE 

F{d)^-smz, z-^-FT- (30) 

The factor -F(O) = 1 and it decreases with distance. This means that the contribution of 
remote structures to the integral (29) is suppressed. The width of the first peak oi F[d) 

la = (31) 

(corresponds to z — tt) determines the attenuation length: at d < la the effects of structures 
are not suppressed. The better the energy resolution, the larger the attenuation length, and 
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consequently, deeper structures can be seen by the neutrino "microscope". This explains, 
e.g., why for the solar neutrinos the zenith angle dependence of the Earth mater effect is fiat 
and there is no enhancement of the regeneration for the core crossing trajectories in spite of 2 
- 3 times larger density. Indeed, for the solar neutrinos with £■ ~ 10 MeV and AE/E = 0.3, 
we obtain = 1000 km and therefore the contribution of the core is attenuated. On the 
contrary, small structures (~ 10 km) near the surface can produce strong effect. The 
attenuation length increases with energy. For ~ 25 MeV and AE/E = 0.2, = 4000 
km, so that the core of the Earth can be probed by the SN neutrinos [42]. 

Another insight into phenomena can be obtained using the adiabatic perturbation theory 
which leads to 



freg = e(i?) sin^ 26'sin^[0.5$„i(xo Xf)] + sm29Re[c{xo Xf)]. (32) 
Here e{R) is the expansion parameter at the surface of the Earth and 

c(x) = r dx'^^^^e'^-^^'-^^^ (33) 
Jxo dx 

is the amplitude of transition between the eigenstatcs in matter (the adiabaticity violation 
effect). In the adiabatic case, c{x) = 0, the second term in (32) is absent. The adiabaticity 
condition is broken at the borders of shells only. Due to sharp density change we have 
for the jth border: d9m{x)/dx\j = AVjS{x — xj), where AVj is the jump of the potential 
between the shells. The integration in c{x) is trivial and simple computations give [20] 

2Esin2 2^ . <I>o V- A.. . ^7 

Here $o and $j are the phases acquired along whole trajectory and on the part of the 
trajectory inside the borders j. This formula corresponds to symmetric profile with respect 
to the center of trajectory. Using (34) one can easily infer the attenuation effect. 



5 Conclusions 

The large mixing MSW conversion provides the solution of the solar neutrino problem: it 
leads to determination of Am^2 ^'^d 612. Now we have detailed physical picture of the 
conversion and its very precise analytical description both inside the Sun and in the Earth. 

Interesting relation emerges between 612 determined in the solar neutrinos and the 
Cabibbo angle: 612 + Oc = 7r/2. If not accidental, it has important imphcation for the 
fundamental physics. 

The small mixing MSW conversion driven by the 1-3 mixing can be realized for the 
supernova neutrinos. Study of these neutrinos will give information on the 1-3 mixing and 
type of mass hierarchy; it opens unique possibility to perform monitoring of shock wave. 
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A number of matter effects can be realized for neutrinos propagating inside the Eartfi: 
(i) the resonance enhancement of oscillations; (ii) the parametric effects in the multi-layer 
medium; (iii) the attenuation effect for the low energy neutrinos. The first two effects can 
be seen in experiments with the atmospheric and accelerator neutrinos. They will play 
important role in determination of the oscillation parameters and establishing the type of 
neutrino mass hierarchy. The attcmiation effect is realized for the solar and supernova 
neutrinos, it describes the loss of sensitivity to remote structures of the density profile. The 
effect is crucial for the oscillation tomography of the Earth. 
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